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Based on experimental results and our previous theoretical work, a microscopic theory of high tem- 
perature superconductivity is conjectured. In this conjecture, superconducting and antiferromag- 
netic long-range orders are driven by interlayer coupling. Strictly in two dimensions, the microscopic 
Hubbard model has an (resonating valence bond) insulator-to-metal transition aX x = Xc near op- 
timal doping for zero temperature, leading to a quantum critical point, and one of the crossover 
lines is given by the pseudogap temperature T* . We argue that various singular and non- Fermi 
liquid properties observed near optimal doping are due to the presence of this quantum critical 
point. In our conjecture, the crossover line T* also practically divides the superconducting region 
into two, depending on the doping level with respect to Xc- For x < Xc the superconducting state 
has significant antiferromagnetic correlations, while for x > it has virtually no antiferromagnetic 
correlations, thus justifying the conventional BCS theory based on the noninteracting electrons. 
Inelastic neutron scattering resonance and systematically reduced superfluid density in the super- 
conducting state below Xc have their natural explanations in the present scenario. The present 
approach supports interlayer pair tunneling model where the superconducting condensation energy 
comes from the lowering of the c-axis kinetic energy in the superconducting state. Comparison of 
the present scenario with some of the leading theories based on the Hubbard and t — J models is 
given. The generic features of both hole-doped and electron-doped cuprates as well as heavy-fermion 
superconductors may be understood in the unified framework within the present picture. 

PACS numbers: 71.10.Fd, 71.27.-ha 
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I. INTRODUCTION 

Since the discovery of high temperature superconduc- 
tivity in copper oxide compounds, Q enormous experi- 
mental and theoretical effort has been made in order to 
understand various anomalous behaviors in normal and 
superconducting (SC) states of these materials. Right 
after its discovery, it was recognized that the high tem- 
perature cuprate superconductors have in common lay- 
ered perovskite-like crystal structures which consist of 
conducting Cu02 planes separated by layers of other el- 
ements. These latter layers function as charge reservoirs 
and mobile charge carriers (holes or electrons) supplied 
from them are believed to reside mainly within the Cu02 
planes. Through extensive worldwide effort as well as the 
improvement of sample quality over the years, many con- 
sensuses have been reached in experimental side. First, 
let us start by showing a generic phase diagram (Fig. 
|l|(a)) of a hole-doped cuprate Ln2_£i;Sr;j,Cu04 in the dop- 
ing {x = 1 — n) and temperature (T) plane. Q Since 
this compound has only one Cu02 layer in a unit cell, 
it is one of the best materials to study physics taking 
place in a two-dimensional Cu02 plane. In other hole- 
doped compounds, the presence of multi-layers and CuO 
chains makes it more complicated to extract the intrinsic 
features associated with a single Cu02 plane. But the 



Cu02 plane in Ln2--2;Sr2;Cu04 is still not totally isolated 
from other Cu02 planes in different unit cells. 

Near half-filling and at low temperature, antiferromag- 
netic (AF) long-range order appears with T/v = 250-300 
K at a; = 0. It is destroyed by 2% doping concentration. 
When X reaches 0.06, SC long-range order starts to ap- 
pear, and it is also destroyed by 30% doping. In between 
them, Tc reaches a maximum value of 40 K at a; ~ 0.16. 
The SC gap was found to have mainly d-wave charac- 
ter with possibility of a small mixture of other angular 
momentum states, in contrast to conventional BCS 
superconductors with an isotropic s-wave gap. The 
low temperature phase between T/v and Tc is often des- 
ignated as spin glass (SG) phase. 

Various recent experiments also show the existence of 
a crossover temperature T* larger than Tc in a doping 
range of x = to a; ~ 0.18 - 0.19. Below this 

pseudogap temperature T*, the low frequency spectral 
weight begins to be strongly suppressed. Surprisingly 
the doping dependences of T* and Tc are completely dif- 
ferent in spite of their close relationship suggested by an- 
gle resolved photoemission (ARPES) iQjs), tunnehng |l^] 
and NMR experiments, jll] At optimal doping where Tc 
is maximum, various non-Fermi liquid (NFL) properties 
are observed in the normal state. These include the lin- 
ear temperature dependence of ab-plane resistivity, the 
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quadratic T dependence of Hall angle and so on up to 
1000 K. Far beyond optimal doping, the normal state 
properties are well described by the conventional Lan- 
dau Fermi liquid. Several anomalous behaviors have been 
observed in the SC state as well. Near optimal doping 
and at underdoping, a sharp resonance (almost energy 
resolution limited) is observed in inelastic neutron scat- 
tering experiments, [ p^ Furthermore in the same doping 
range, the superfluid density is systematically sup- 
pressed with decreasing doping in spite of increasing SC 
gap amplitude. In the overdoped regime, however, the 
SC properties appear to be well explained by the con- 
ventional weak coupling BCS theory. 

Although they are not usually addressed in the con- 
text of a phase diagram, the following experiments de- 
serve special attention in order to understand the com- 
plete picture of the physics. A recent ARPES experi- 
ment for an insulating cuprate Sr2Cu02Cl2 [0 showed 
nearly isotropic and quite similar band dispersions along 
(7r/2, 7r/2) - (tt, 0) and (7r/2, 7r/2) - (0, 0) directions. Fur- 
thermore a c?-wave-like modulation of the insulating gap 
in Ca2Cu02Cl2 jl^ requires even at half-filling the pres- 
ence of strong pairing fluctuations. In another important 
experiment, a 61 Tesla pulsed magnetic field suppressed 
superconductivity in Ln2-a;Sr2;Cu04 single crystals and 
revealed an insulator-to-metal (IM) crossover for both 
ab-plane resistivity pab and c-axis resistivity pc near opti- 
mal doping, [ p^ The latter experiment may shed some in- 
sight into the underlying physics of normal state cuprates 
which is usually masked by the presence of SC long-range 
order. Any successful theory should explain all these fea- 
tures in a natural and yet unified way. 

As to the theoretical side, right after the discovery of 
high temperature superconductors, Anderson first 
proposed the one-band Hubbard model as the simplest 
Hamiltonian which might capture the correct low energy 
physics of copper oxides. Zhang and Rice jl^ also de- 
rived an effective t — J Hamiltonian from the more re- 
alistic three-band Hubbard model. The t — J Hamilto- 
nian was already known to be the large U limit of the 
Hubbard Hamiltonian under certain assumptions. The 
Hubbard model is described by the Hamiltonian in which 
Ci^CT destroys an electron at site i with spin u on a two- 
dimensional square lattice 

H ^-t J2 4<tCj\'^ + ■ (1) 

i is a hopping parameter between nearest neighbors («, j) 
and U denotes local Coulomb repulsion. It is believed 
that a realistic strength of the Coulomb repulsion lies in 
between the weak and strong coupling regimes, namely, 
U W — 2W where W is the bandwidth of 8t in two 
dimensions. So far many theories have been proposed 
to understand the various anomalous properties of the 
cuprates by directly invoking the Hubbard model or its 



variants (one-band and three-band Hubbard, and t — J 
models), or other phenomenological models. However, at 
present there is no consensus on which kind of approach 
is most appropriate for a given Hamiltonian or even on 
which Hamiltonian is most relevant for the cuprates. At 
the end of this paper, comparison of our approach with 
some leading theories based on the Hubbard and t — J 
models will be given. 

II. REVIEW OF OUR PREVIOUS WORK 

First, let us begin by summarizing the main results of 
our previous work and presenting a proposed phase 
diagram in order to set the stage for the present study. 
As a first step to the microscopic theory of high-Tc super- 
conductivity, it is of great importance to identify pairing 
interaction which may eventually lead to a SC instabil- 
ity at low temperature, just like the Cooper problem in 
the development of the BCS theory. In a recent study, 
we examined the possibility of extracting pairing inter- 
action directly from the Coulomb repulsion itself without 
the exchange of bosonic degrees of freedom such as spin 
fluctuations. 

We found that among several gap symmetries pairing 
interaction with the c?a:2-j;2 {(f)d{k) = cos fca, — cos fcj, ) sym- 
metry is most strongly induced from the local Coulomb 
repulsion U . It is consistent with a weak coupling renor- 
malization group (RG) calculation. Due to its direct 
Coulombic origin, this mechanism for pairing correlations 
is similar to the lattice version of the Kohn-Luttinger the- 
orem |^l[]. However, the details are somewhat different. 
The Kohn-Luttinger theory requires a sharp Fermi sur- 
face and a resulting long-range oscillatory behavior in 
real space of the effective interaction, while our approach 
does not. 

AF correlations are also induced from the local 
Coulomb repulsion at the same time as d-wave pairing 
correlations are. Both the induced pairing and AF cor- 
relations are found to increase with decreasing doping. 
The simultaneous induction of pairing and AF correla- 
tions from the same Coulomb repulsion may be under- 
stood in the context of the t — J model. Projecting out 
the doubly occupied sites in the large U limit of the Hub- 
bard model generates or induces the Heisenberg term in 
the t — J model. It has in general both pairing and AF 
correlations. The increasing pairing and AF correlations 
with decreasing doping found in our previous work is also 
realized in the t ~ J model (in a relative sense) as de- 
creasing hopping probability with decreasing doping. In 
that paper, we argued that the induced pairing correla- 
tions are the microscopic origin of pseudogap behavior in 
underdoped cuprates. This is because the pairing fluctu- 
ation origin of the pseudogap and its doping dependence 
are consistent with many experiments. This is also be- 
cause of the excellent agreement of our calculations with 
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the observed energy dispersion of the insulating cuprates, 
SraCuOzClz |l| and CasCuOaCla. [|| It is expected 
that the high energy pseudogap at half-fiUing continu- 
ously evolves into a relatively smaller pseudogap away 
from half-filling. 

Based on the doping dependences of the induced pair- 
ing and AF correlations, a phase diagram was pro- 
posed for the high temperature superconductors, which is 
shown in Fig. |l|(b). First, far away from half- filling the 
short-range correlations due to local Coulomb repulsion 
are very weak and don't produce strong enough pairing 
interaction so that Tc vanishes. With decreasing doping, 
they start to induce both pairing and AF fluctuations, 
but the latter is probably not strong enough to compete 
with the former. Thus up to optimal doping level, local 
Coulomb repulsion plays as a driving force to supercon- 
ductivity and Tc keeps growing. When electrons enter 
into the underdoping regime, the induced AF fluctua- 
tions start to play as a SC phase coherence-breaker by 
creating locally the spin density wave (SDW) state. It 
strongly breaks time-reversal symmetry and causes Tc to 
decrease. As far as T* is concerned, beyond optimal dop- 
ing, induced AF fluctuations are weak and thus the pseu- 
dogap temperature T* closely follows Tc. The strength 
of induced pairing interaction increases with decreasing 
doping. This makes T* increase, creating a huge pseudo- 
gap region between T* and Tc. 

Near half-filling and at low temperature, AF long- 
range order is stabilized by the three-dimensional effect. 
The AF spin fluctuations become strong more or less 
around the AF phase boundary, leading to another kind 
of a crossover temperature Trvb. Below Trvb both 
AF and pairing fluctuations are strong and coexist. It is 
expected that the SC long-range order completely van- 
ishes, when it enters into the RVB-like region. The T = 
ground state between the AF and SC phase boundaries 
(in a clean sample) is the RVB insulator in which an in- 
sulating gap is identified with the spingap A. Note that 
Tc, T* , Tn and T^yb are determined as a result of close 
competition between the induced pairing and AF corre- 
lations and that any imbalance of the two correlations 
may drastically change the phase diagram. The pseu- 
dogap behavior in this scenario is as generic as AF cor- 
relations, because they are simultaneously induced from 
the same local Coulomb repulsion. In our picture T* is 
different from a temperature below which the critical be- 
havior (0(2) symmetry) of SC long-range order is 
found. In this previous work, however, interlayer cou- 
pling effect and the possibility of existence of a quantum 
critical point (QCP) were not considered. 



III. RVB INSULATOR-TO-METAL TRANSITION 
NEAR OPTIMAL DOPING AT ZERO 
TEMPERATURE: QCP IN THE 
TWO-DIMENSIONAL HUBBARD MODEL 

In principle, the above phase diagram resulting from 
the close competition between the induced pairing and 
AF correlations should be quantitatively similar for differ- 
ent cuprates with similar C/'s. This is because these corre- 
lations are completely determined by the same Coulomb 
repulsion U at given x and T . For most of the hole- 
doped cuprates, J = At'^/U lies in the range of 120-140 
meV, indicating quite similar C/'s with 20% variation at 
most. One the other hand, Tc's, which are expected to 
be similar with 20% difference, widely vary from 40 K for 
Ln2-2;Sr2;Cu04 to 133 K for IIgBa2Ca2Cu308. This par- 
ticular feature strongly suggests that the true SC long- 
range order is driven by interlayer coupling. It is consis- 
tent with the general trend that the cuprate compounds 
with more Cu02 planes in a unit cell show higher Tc. In 
the present paper, the interlayer coupling means not only 
the coupling between Cu02 planes in a unit cell, but also 
between other Cu02 planes in different unit cells. Strictly 
in two dimensions, AF long-range order cannot appear at 
finite temperature due to the Mermin- Wagner theorem 
, while SC long-range order can be stabilized at finite 
temperature by the Kosterlitz-Thouless phase transition. 

n 

However, there is possibility that SC long-range or- 
der is also destroyed at finite temperature mainly by the 
phase coherence-breaking effect of the induced AF fluctu- 
ations and partly by quantum fluctuations. This is con- 
sistent with the rigorous result by Su and Suzuki, p^ ] 
These authors proved the nonexistence of dx2-y2 super- 
conductivity (long-range order) at any nonzero temper- 
ature in the two-dimensional Hubbard model. Note that 
the proof by Su and Suzuki does not exclude increasing 
pairing correlations [ pTf with decreasing temperature just 
like the AF correlations for a half-filled Hubbard band. 
There is also possibility that strictly in two dimensions 
both SC and AF long range orders are absent even at 
zero temperature. This is because AF fluctuatings with 
spin triplet character (pairing fluctuatings with spin sin- 
glet) destroy SC long-range coherence with spin singlet 
character (AF long-range coherence with spin triplet). In 
a sense, the induced pairing and AF correlations may act 
in two dimensions like frustrating the long-range coher- 
ence of the other correlations. Unlike the SC long-range 
order, the absence of AF long-range order at zero tem- 
perature is not so crucial for our argument. Thus in our 
first conjecture, only by the three-dimensional effect are 
stabilized the AF and SC long-range orders at finite tem- 
perature and even at zero temperature. 

Now in order to understand the intrinsic electronic 
properties in a totally isolated Cu02 plane, let us do a 
thought experiment on the phase diagram (Fig. |l|(a)) by 
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slowly turning off the interlayer coupling. It is clear that 
slow turning off of the interlayer coupling reduces the AF 
and SC phase boundaries. At the same time the region 
of the RVB insulating ground state expands, as shown in 
Fig. |2|(a). When the interlayer couphng is completely 
disconnected from neighboring Cu02 planes in the same 
unit cell or in other unit cells, the AF and SC long-range 
orders totally disappear. And the RVB insulating ground 
state persists from half-filling all the way up to x = Xc 
near optimal doping, as shown in Fig. |^(b). Consid- 
ering T* meets Tc at slightly overdoping region and its 
curvature, x = Xc should be somewhat larger than opti- 
mal doping (xc — 0.18 — 0.19). Since T* involves a short 
distance scale or rather a high energy scale, T* obtained 
in two dimensions is not expected to change substantially 
by a small interlayer coupling. 

Thus in our second conjecture, the two-dimensional 
Hubbard model has an RVB insulator-to-metal transi- 
tion at X = Xc for T = 0, leading to a QCP, and one of 
the crossover lines is given by the pseudogap tempera- 
ture T*. This feature is presumably what Boebinger et 
al. ]l6t found in their experiment of an insulator-to-metal 
crossover near optimal doping. As is well known from an 
extensive study of QCP by Sachdev, an immediate 
consequence from the presence of QCP is to create two 
crossover lines which occur at T ~ lAj^"^, where z is the 
dynamical exponent and v is the correlation length ex- 
ponent (Fig. 1(c)). The crossover line above the RVB 
insulating ground state is identified with the pseudogap 
temperature T*. The state above the QCP, sandwiched 
by the two crossover lines, is described by placing the 
A = scale-invariant critical theory at nonzero temper- 
ature. This is probably the region where various singular 
and NFL properties have been observed in the nor- 
mal state and phenomenological marginal Fermi liquid 
(MFL) theory by Varma et al. works well. The fact 
that the pseudogap temperature T* naturally divides the 
pseudogap region and the NFL region, justifies a com- 
monly used experimental method [|3l| of determining the 
pseudogap temperature. For instance, in dc resistivity 
measurements, T* is a temperature where the resistivity 
starts to deviate from a linear temperature dependence. 

Since we have identified strictly in two dimensions the 
QCP of the RVB insulator-to-nietal and concomitant 
NFL behaviors in optimally doped cuprates, let us turn 
on the interlayer coupling in order to find the realistic 
phase diagram in quasi- two dimensions. Before doing 
that, it is important to point out the differences in elec- 
tronic properties below and above Xc at low temperature. 
In our third conjecture, at low temperature (T* ~ 0), the 
crossover line T* roughly divides the doping-temperature 
plane into two regions, one with AF correlations and the 
other without them. The identification of T* at low tem- 
perature as a crossover of the AF correlations may be 
understood by the following two observations. 

In our previous work, |19| another crossover line Tjiy b , 



which is not specifically addressed in this paper, was de- 
fined as a line where the AF correlations are just strong 
enough to destroy the SC long-range phase coherence, as 
shown in Fig. |l|(b). When the interlayer coupling is 
completely turned off, it is easy to see Tbvb approach 
T* at low temperature. As another supporting evi- 
dence, we and other groups |3^j3^ found that most of 
the low temperature region including the SC state below 
X ~ 0.2 is inside of the AF phase in a mean-field study 
of the t — J model. In a mean-field approximation, long- 
range order already sets in when the corresponding cor- 
relation length reaches roughly one lattice spacing. This 
dictates the above mean-field result to be interpreted as 
the presence of short-range AF correlations for x < 0.2 
at low temperature. If this conjecture turns out to be 
correct, there is possibility that at low temperature and 
near a; ~ Xc a band Kondo effect is realized in finite di- 
mensions, as is done in infinite dimensions. |^6| In the 
band Kondo effect, a nonvanishing local magnetic mo- 
ment at site j, {Sf), due to local Coulomb repulsion is 
effectively screened out by surrounding (but the same 
kind of electrons at site i) itinerant electrons. 

From a purely pairing correlation point of view, the 
pairing correlations become singular with decreasing 
temperature at any doping concentration. This makes 
the crossover region of the pairing correlations much 
wider than that of the AF correlations at low temper- 
ature near x ~ Xc. As another supporting evidence on 
this point, we again invoke our previous mean- field phase 
diagram. The mean-field Tc line, which is interpreted 
as a crossover of short-range pairing correlations, extends 
up to X ~ 0.35 beyond x = Xc where the mean-field AF 
order vanishes. It makes the crossover region of pairing 
correlations broad with respect to x = Xc at low temper- 
ature. 

In this situation, a potential location with the high- 
est Tc, when the interlayer coupling turned on, is near 
X = Xc in which the phase coherence-breaking AF corre- 
lations nearly vanish but the pairing correlations are still 
robust. The resulting phase diagram (Fig. 0(a)) will 
look like one where T* falls from a high value onto the 
Tc line rather than the other where T* smoothly merges 
with Tc in the slightly overdoped region, as recently ar- 
gued by Tallon and Loram. p7| In our scenario Tc is 
never part of the T* line. For different compounds, dif- 
ferent strength of interlayer pairing hopping drives SC 
long-ranger order in the same background of the two- 
dimensional electron system. This leads to a universal 
relation |3|] Tc/T^^^ = 1 - 82.6(x - 0.16)^, when Tc 
is scaled by T^^-^. Away from half-filling, the AF cor- 
relations manifest their existence most strongly in the 
SC state, because they easily destroy the SC long-range 
phase coherence. Our scenario also predicts that due to 
scatterings with AF fluctuations, quasiparticle scattering 
rate remains finite for x < Xc even in a clean sample and 
at T = and it is universal (up to the variation of J or 
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U for different compounds). For a; < a;^, in our picture, 
the Landau quasiparticle with nonvanishing quasiparticle 
residue does not exist in the normal state due to strong 
scatterings with pairing and AF fluctuations or strong in- 
fluence from the QCP. But it can be stabilized in the SC 
state (T < Tc) where its coherence is restored through 
a pair hopping process along the c-axis. Because of the 
strong influence from the QCP, the NFL takes a special 
structure near optimal doping, presumably the marginal 
Fermi liquid by Varma et al. [pO[ 



IV. INELASTIC NEUTRON SCATTERING 
RESONANCE AND REDUCED SUPERFLUID 
DENSITY IN THE SC STATE 

As long as the strength of pairing correlations is con- 
cerned, it lies in the intermediate to weak coupling 
regimes near optimal doping, as shown in Ref. |19|. If a 
screening effect is correctly taken into account in obtain- 
ing Vind (by considering particle-particle and particle- 
hole channels on equal footing), Vind can be further re- 
duced near and beyond optimal doping. Thus the cru- 
cial point in proper understanding of the SC state is not 
the validity of the weak coupling BCS theory, but more 
importantly whether Cooper pairs are constructed from 
antiferromagnetically correlated electrons or not. In fact 
depending on the doping concentration with respect to 
a; = Xc, the SC state can be qualitatively different. For 
X < Xc the SC state has significant AF correlations, while 
for a; > a;c it has virtually no AF correlations, thus jus- 
tifying the conventional BCS theory based on the nonin- 
ter acting electrons. 

First, we discuss a strong inelastic neutron scatter- 
ing resonance Q observed in the SC state of the un- 
derdoped and optimally doped cuprates (a; < Xc)- The 
striking feature of a sharp (almost energy resolution lim- 
ited) resonance peak is that its energy E,. does not 
shift towards lower energy when approaching Tc, but its 
intensity vanishes upon heating at T^. On experimen- 
tal grounds, the resonance peak is often considered as a 
collective mode in the SC state. The inelastic neutron 
scattering experiment measures the dynamical spin-spin 
correlation function 



S(q,uj) 
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where |0) and \n) are the ground and excited states of the 
system, and a;„o is the energy difference between these 
states, and St = Y^rci ^ cr , . 

In the present scenario, the neutron scattering reso- 
nance is caused by the presence of AF correlations in 
the SC state for a; < a;c. In a recent mean- field study 
of the interplay between antiferromagnetism and d-wave 
superconductivity, pO[ we noted that when d-wave SC or- 
der and AF order coexist another order parameter 



with a spin-triplet, (c^:^^ -^c ^: ^), is dynamically gener- 
ated. Q is the AF wave vector (tt, tt) in two dimensions. 
This is a generic feature for fermionic systems and can- 
not be obtained in a purely bosonic description of pairing 
and AF correlations. As long as the two correlations are 
strong (but both of them need not be in long-range or- 
dered states), the spin-triplet correlations can be robust. 
Since both pairing and AF correlations are present below 
Xc as noted in the previous Section, in general the ground 
state has three different (pairing, AF and spin-triplet) 
correlations. To have a maximum transition amplitude 
(n 1 1 0)1, one finds two possible ways of transition 
between the three correlations, namely, AF ^ AF and 
pairing ^ spin-triplet. The former is most effective near 
half-filling, due to the presence of the AF long-range or- 
der in that region. On the other hand, the latter is most 
operative near x = Xc, because the effective strength of 
SC long-range order is strongest there. |0 Between at 
half- filling and at a; = Xc, the two ways of transition 
compete. 

Near a; = a;c where a strong neutron resonance was ob- 
served and the second way of transition is most effective, 
the resonance energy is given by twice of the maximum 
gap energy, and its intensity by the product of d-wave 
SC order and SDW order parameters Since there 

are only short-range AF correlations in the SC state, the 
resonance intensity is roughly the product of superfluid 
density and effective strength of AF correlations. In op- 
timally doped and underdoped samples {x < Xc)^ the SC 
gap amplitude is already fixed as the pseudogap size be- 
low T* and thus the resonance energy does not shrink 
towards lower energy when approaching Tc- However, its 
intensity vanishes (or becomes strongly suppressed) when 
temperature is increased toward Tc, following an order 
parameter-like behavior due to its proportionality to the 
superfluid density. This suggests that the sharp inelastic 
neutron resonance near optimal doping is a direct con- 
sequence of the collective excitations of the spin-triplet 
ground state. Observed incommensurate magnetic peaks 
may come from a band structure effect at the Fermi en- 
ergy. There is also possibility that they are from the for- 
mation of inhomogeneous stripe structure. p2[ It is im- 
portant to note that the (energy resolution limited) col- 
lective mode nature of the neutron scattering resonance 
is best understood, when the pairing and AF correlations 
are treated on equal footing. 

Recently Demler and Zhang have reached a similar 
result to ours in the context of SO (5) symmetry, although 
details are somewhat different. In SO (5) theory, the spin- 
triplet amplitude is not an order parameter, but the gen- 
erator of infinitesimal rotations between AF and SC or- 
der parameters. In the mixed AF-f SC-|-spin-triplet phase 
that is discussed here, the spin-triplet state acquires a 
nonzero order parameter in the ground state. The above 
argument which is valid near a: = a;c drastically changes 
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with decreasing doping. In underdoping region where the 
first way of transition starts to become more operative, 
the resonance energy is no longer given by twice of the 
maximum gap energy. Otherwise, it would keep increas- 
ing with more underdoping, because the pseudogap size 
increases with decreasing doping. The intensity also does 
not behave like an order parameter any more, as it does 
near optimal doping. As the doping is decreased, the 
resonance energy becomes soft and becomes a Goldstone 
mode of the AF order near half-filling. To our opinion, 
there is no obvious reason for linear scaling of the reso- 
nance energy with Tc- At present it is uncertain whether 
this linear scaling is just a coincidence or comes from a 
deeper theoretical origin. It is expected that at optimal 
doping region or in slightly underdoping region two res- 
onance peaks can appear in inelastic neutron scattering 
experiments, one from the first way of transition and the 
other from the second way. 

With decreasing doping, the induced pairing correla- 
tions as well as the phase coherence-breaking AF cor- 
relations increase in underdoped and optimally doped 
samples. As a result, with decreasing doping, the SC 
gap amplitude increases but superfluid density or Tc de- 
creases. The resulting { Ad)ma.x/ KbTc ratio is strongly 
doping dependent, monotonically increasing with de- 
creasing doping below Xc- Above Xc (overdoping), how- 
ever, it is expected that the ratio approaches more or 
less the BCS mean-field value. This feature cannot be 
understood in the absence of AF correlations which are 
allowed in the model. The effective strength for the SC 
long-range order is also strongly doping dependent and is 
largest near Xc- It decreases below Xc due to the increas- 
ing phase coherence-breaking AF correlations and also 
above Xc owing to the decreasing pairing correlations, 
as discussed in the previous Section. In this respect, it 
is not surprising to find that the superfluid density and 
the SC condensation energy have their maximum values 
near Xc, and decrease below and above Xc- In the 
present scenario, the pseudogap is virtually unchanged 
by an applied magnetic field, because the characteristic 
energy scale for the pseudogap. A, is much larger than 
the Zeeman energy. On the other hand, the SC long- 
range order is relatively easily destroyed by it because of 
its phase coherence-breaking nature. 

Finally it is worthwhile to comment on the SC conden- 
sation energy. In the present scenario the SC long-range 
order is stabilized only through a pair hopping process 
along the c-axis (due to the interlayer coupling). It forces 
the SC condensation energy to come from the lowering 
of the c-axis kinetic energy in the SC state. This inter- 
layer coupling theory was already proposed by Anderson 
and others, and many features are consistent with c- 
axis optical measurements. p5| Recently some significant 
discrepancies found in the measured c-axis penetration 
depth in Tl2Ba2Cu06+5 and in the prediction from the 
interlayer tunneling model, were resolved by Chakravarty 



et al. |46| by subtracting fluctuation effects in the elec- 
tronic specific heat data. 



V. GAP SYMMETRY IN THE SC STATE 

In this Section, we discuss the gap symmetry of the SC 
order parameter for cuprates. In the previous study, p9| ] 
for several d-wave type symmetries the effective strength 
of pairing correlations, (A^(O)Ag(O)), was found as 



AN 



Y A(fc) 



2UN 



^^sdw and 



where X{k) = ^ {{e{k) - e{k + Q)) /2f 

E±{k) = (e(fc) + e{k + Q))/2 ± A(fc). Generally Ag(i) 
is defined as Ag(i) = \ Y.5 9{5){ci+SACi^[ - Ci+5,j,Ci,|), 
where g{6) is an appropriate gap structure factor in real 
space. e(fc) = —2t{coskx + cosfcy) — ^ for nearest neigh- 
bor hopping, /i is the chemical potential controlling the 
particle density n, N the total number of lattice sites, 
f{E) the Fermi-Dirac distribution function, 4>g{k) the 
Fourier transform of g{5)^ and the summation accompa- 
nied by the prime symbol is over wave vectors in half 
of the first Brillouin zone. For local s-wave and ex- 
tended s-wave symmetries, there are additional contri- 
butions to the above equation. But the second term in 
Eq. ^ is still the major factor determining whether or 
not (A+(0)A,(0)) > (A+(0)A,(0))o. 

Purely from a symmetry reason are ruled out sev- 
eral gap symmetries such as local s-wave {4>{k) = 1) 
and dxy {^'{k) = 2 sin fc^; sin /cj,), both of which have 
4>{k -\- Q) = 4'{k). In the real space representation, a 
pair with a local s-wave symmetry involves an upspin 
(downspin) electron and a downspin (upspin) electron at 
the same site i. A pair with a dxy symmetry the combi- 
nation of an upspin (downspin) electron in site i and the 
downspin (upspin) electrons in the next-nearest neigh- 
bors with an appropriate sign. Thus the pair configura- 
tions of the above two symmetries are directly against 
that dictated by the strong local Coulomb repulsion U . 
Eventually the most stable gap symmetry among various 
possible types with <l)[k + Q) = — 0(fc), is determined by 
the band structure near the Fermi surface. Near half- 
filling |(/'(A:)| of dx2^y2 symmetry {4>{k) = coskx — cosfcy) 
is much larger than, for instance, that of extended s-wave 
{4>{k) = coskx + cos ky) near the Fermi surface. It en- 
ables dx2-y2 gap symmetry to be realized in most cases. 
Far away from half-filling, however, the extended s-wave 
form factor becomes important and there is possibility 
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that the SC gap contains a substantial fraction of ex- 
tended s-wave component, leading to a dx'^-yi + is gap 
symmetry. 

Recently there is some controversy over the gap 
symmetry in electron-doped cuprates Nd2-a;Cej;Cu04_5 
and Pr2_a;Ce2:Cu04_5, namely, isotropic s-wave against 
anisotropic c?-wave. | |47| According to the above princi- 
ple of finding the most stable gap symmetry, local (or 
isotropic in momentum space) s-wave is ruled out from 
the beginning. It is due to its incompatibility with the 
underlying strong local Coulomb repulsion. The correct 
gap symmetry for electron doped cuprates should be ex- 
tended s-wave or d^^-y^ or a mixture of these symme- 
tries, depending on the projected density of states at the 
Fermi energy for a given symmetry. The gap of an ex- 
tended s-wave symmetry vanishes only at the AF zone 
boundary. According to ARPES experiments |Q and 
band structure calculations, [^9| however, the Fermi sur- 
face of electron doped cuprates is likely to be circular 
near the zone center. In this situation a gap node of 
the extended s-wave symmetry does not appear in the 
Fermi surface. Consequently this symmetry, if it is cor- 
rect, might have been interpreted as the isotropic s-wave 
in previous experiments. This feature may imply that 
the observed pairing symmetry itself does not tell us any 
crucial information about the mechanism of supercon- 
ductivity. 



VI. COMPARISON WITH SOME LEADING 
THEORIES BASED ON THE HUBBARD AND 
T J MODELS 

In this Section we compare the present approach with 
some leading theories for the high-T^ superconductiv- 
ity. Perhaps the closest theories to our approach are 
those based on the Anderson's RVB state in one way 
or another. Anderson and his co-workers first applied a 
mean-field approximation to the t — J model. The 
mean-field theory studied by these authors and by oth- 
ers |5l|-|5^] has been a starting point for further devel- 
opment of the theory such as 1 /N expansion theory |5^] 
and gauge theory |^,^ of the t — J model. In slave 
boson theory of the t — J model, typically two mean-field 
order parameters are considered 



(4) 



together with (bi) In the slave boson representation, a 
physical electron is decomposed into a spinon (ferniion) 
and a holon (boson), c/"^ = f^^bi. 

Depending on the vanishing or nonvanishing of Ay and 
{bi), the doping and temperature plane is divided into 
four regions. Region I with Ay = and {bi) 7^ is a 
Fermi liquid phase. Region II with A^ ^ and {bi) — 



is the spingap phase, in which a d-wave gap appears in the 
fermion spectrum without Bose condensation of holons. 
Region III with A^- 7^ and {bi) ^ indicates SC long- 
range order in physical electrons. Region IV with Ay = 
and {bi) = is designated as the strange metal phase, 
because it shows various non- Fermi liquid features. 

In many respects, the slave boson mean- field theory of 
the t — J model [^-56| has shed some important insight 
into the microscopic understanding of the cuprate super- 
conductors. This is because the predicted phase diagram 
is, at least, qualitatively consistent with experiments, and 
the pseudogap is closely related to a spingap, and further- 
more it starts from the microscopic model as opposed to 
other phenomenological models. However, there are also 
some serious problems with the slave boson mean-field 
theory, as noted by Ubbens and Lee. |5|] 

One of them is that the temperature scale for Bose con- 
densation of holons is too high. Furthermore the maxi- 
mum Tc, which is determined by the two lines A^ ^ 
and {bi) / 0, occurs at too small doping concentration 
{x < 0.06). At this doping level, the SC long-range order 
even does not appear in cuprate superconductors. Close 
to half-filling, several exotic phases have been reported to 
be stable such as mixed phases ]52| ( equivalently 7r-flux 
phases p3|) , dimerized phases, [ ^3| , MI ^^'^ staggered flux 
phases. [65 56 1 It is unclear whether these states are real- 



ized or not in cuprates. In a recent paper ]33| we argued 
that these problems can be naturally resolved, when AF 
correlations, the weaknesses of a mean-field approxima- 
tion, and the limitation of the t — J model near half- filling 
are properly taken into account. 

In the original RVB theory, Anderson sug- 

gested that the RVB ground state may be obtained by 
Gutzwiller projection on the BCS ground state. As 
shown in the previous study, |ll9| however, strong local 
Coulomb repulsion U induces the AF and pairing correla- 
tions at the same time. This makes one Gutzwiller pro- 
jection on the noninteracting electrons enough to yield 
the RVB ground state. Unless the interlayer coupling is 
introduced in the RVB based theories, they suffer from 
having only one energy scale A which has to explain T* 
and Tc at the same time. Some studies on the t ^ J 
model at half-filling or the Heisenberg model showed that 
a mixed phase [s^ ] (equivalently a tt-Aux phases ||53|| ) is 
the ground state of the model. However, the observed 
energy dispersion in the insulating cuprates Sr2Cu02Cl2 
and Ca2Cu02Cl2, is different from what the 
flux phase predicts. The observed band dispersions along 
(7r/2,7r/2) - (7r,0) and (7r/2,7r/2) - (0,0) are not iden- 
tical and the energy dispersion appears quadratic near 
(7r/2,7r/2) point instead of linear. In fact the experi- 
ments are more consistent with our calculations based 
on the Hubbard Hamiltonian . Although these differ- 
ences (the Hubbard vs. t — J Hamiltonians) are rather 
quantitative, it is not clear whether in some subtle issues 
associated with our conjectures these Hamiltonians give 
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the same answer or not. 

The mechanism to superconductivity proposed in our 
previous and this papers is fundamentally different from 
the conventional one in which superconductivity is driven 
by exchange of some bosonic degrees of freedom such as 
phonons, spin and charge fluctuations and so on. It sug- 
gests that the exchange of spin fluctuations cannot 
be the leading mechanism to high temperature super- 
conductivity, at least, in the Hubbard type model. This 
is because this approach completely neglects the lead- 
ing pairing interaction directly induced from the local 
Coulomb repulsion. But instead it depends on the resid- 
ual pairing interaction from the AF correlations which 
were already induced (simultaneously with the pairing 
correlations) from the Coulomb repulsion. 

VII. NUMERICAL STUDIES FOR THE 
HUBBARD AND T J MODELS 

There have been several numerical studies on the Hub- 
bard and t — J models in two dimensions. []6l| Among 
many important issues, the existence of pairing corre- 
lation which may lead to a SC instability at low tem- 
perature is of particular importance. Several methods 
including exact diagonalization (ED), quantum Monte 
Carlo (QMC) and density matrix rcnormalization group 
(DMRG) support with decreasing temperature growing 
d-wave pairing correlations away from half-filling. This 
is consistent with a weak coupling RG study of the Hub- 
bard model. pO| We believe what many numerical meth- 
ods such as ED, QMC and DMRG find is the pairing 
correlations induced from the local Coulomb repulsion in 
the presence of the induced AF correlations. However, 
the fact that in the numerical calculations the pairing 
strength increases with decreasing doping just as the AF 
correlations do, might have given a misleading impression 
that the pairing correlations come from the exchange of 
spin fluctuations. 

On the other hand, a constrained path quantum Monte 
Carlo (CPQMC) method j6|] appears to show a neg- 
ative answer to the existence of pairing correlations in 
the Hubbard model. This result is in direct conflict 
with that of the RG study. Although it is valid in 
the weak coupling limit, the RG result cannot be qual- 
itatively wrong. It is because of its systematic nature 
and the general experience that many weakly interacting 
systems are more or less smoothly connected to corre- 
sponding strong coupling systems. Thus in the CPQMC 
method, there is possibility that away from half-filling 
their ill-suited trial wave function (SDW solution) gener- 
ates mainly the states with AF correlations but neglects 
those with pairing correlations. Or in the process of con- 
straining paths to avoid the fermion sign problem, the 
d-wave pairing correlations are strongly suppressed. In- 
deed Zhang's calculations may suggest the fate of those 



theories based on purely magnetic correlations such as 
the exchange of spin fluctuations, magnetic polarons and 
so on. 

VIII. HOLE-DOPED CUPRATES AND 
HEAVY-FERMION SUPERCONDUCTORS 

The present scenario may also apply without any dras- 
tic modification to electron-doped high- Tc superconduc- 
tors as well as heavy-fermion superconductors. The main 
qualitative difference between hole-doped cuprates and 
the above compounds ||6^ is that for the latter local mag- 
netic moments from Ce or U atoms residing between con- 
ducting planes play an important role in determining the 
phase diagram and electronic properties. Electrons in 
the Cu02 plane interact with local magnetic moments 
from Ce or U atoms. This enables the AF correlations 
in the Cu02 plane to easily induce AF correlations in 
those atoms, enhancing the overall AF fluctuations in the 
Cu02 plane. Then the close balance between the pairing 
and AF correlations is broken in favor of the latter in the 
conducting electrons. As a result, AF long-range phase 
significantly increases, and accordingly SC phase shrinks 
and remains at the edge of the AF phase boundary. Be- 
cause of the resulting imbalance of the two correlations in 
the conducting plane, there is possibility that the pseu- 
dogap state is completely destroyed or remains at most 
in a narrow region of the x — T plane. The corresponding 
SC ^ap symmetry is determined by the relative sign of 
4>g{k+Q) vs. 4>g{k) and the shape of the Fermi surface. In 
general it can be other than a d-wave symmetry. Thus, 
high temperature superconductors for both hole-doped 
and electron-doped cuprates as well as heavy-fermion su- 
perconductors may be understood in the unified frame- 
work within the present scenario. We also argue that 
differences observed in hole-doped and electron-doped 
cuprates do not come from the different signs of charge 
carriers. They are from the presence of magnetically ac- 
tive atoms (Ce) between Cu02 planes in the latter. If this 
turns out to be correct, the correct model Hamiltonian 
for the heavy-fermion superconductors would be the pe- 
riodic Anderson model or the Kondo lattice model with 
strongly correlated conducting electrons instead of with 
the noninteracting electrons. Various studies of the su- 
perconductivity on the basis of spin fluctuations in these 
compounds and even in organic superconductors should 
be reconsidered for their validity. 

IX. MAXIMIZING THE SC CRITICAL 

TEMPERATURE Tc 

It is also interesting to discuss the issue of maximizing 
the SC critical temperature Tc on the basis of the present 
scenario. There are roughly three ways of increasing T^. 
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First, Tc can be increased by a stronger interlayer cou- 
pling. In the family of Bi2Sr2Ca„_iCu„02n+4 (BSCCO), 
TlBa2Ca„_iCu„02„+3 

(TIBCCO), Tl2Ba2Ca„_iCu„02„+4 (T2BCC0), and 
HgBa2Ca„_iCu„02„+2 (HBCCO) where n =1, 2, 3 and 
4, Tc for larger n is higher than that of smaller n. Except 
for TIBCCO, however, all the others show their highest 
Tc's at n = 3. It indicates the saturation of the interlayer 
coupling strength with increasing the number of layers in 
a unit cell. Making the interlayer region more metallic is 
also suggestive to enhance the interlayer coupling. Sec- 
ond, Tc can be increased by decreasing the local Coulomb 
repulsion U . When HBCCO compound is subjected to 
a high pressure of 30 GPa, Tc increases from 133 K to 
164 K. Applying pressure causes the increase of hopping 
integral between nearest neighbors, effectively decreasing 
the local Coulomb repulsion U . Since Coulomb repulsion 
is also a driving force to the pairing correlations, U can- 
not be arbitrarily small to achieve the highest Tc- There 
must be an optimal strength U (presumably U ~ W) 
in which Tc is maximized. Third, Tc can be also in- 
creased by breaking the close balance between the pairing 
and AF correlations in favor of the former. In principle 
this can be achieved by putting some atoms between the 
Cu02 planes, which respond diamagnetically to the phase 
coherence-breaking AF fluctuations in such a way that 
the overall strength of AF correlations are reduced in the 
Cu02 planes. It is exactly opposite to what happens in 
electron-doped cuprates. Although the implementation 
of this possibility is uncertain in a real situation, this 
last way can be most effective to maximize Tc- 



X. CONCLUSION 

Before closing we comment on a few points. Since the 
conclusion of the present paper crucially depends on the 
three conjectures, it is important to prove or disprove 
them, of course, on the basis of fully systematic or rigor- 
ous calculations. It is believed that the mechanism to su- 
perconductivity and the interplay between antifcrromag- 
netism and superconductivity expounded in this paper 
can be applicable to heavy-fermion superconductors as 
well as organic superconductors. A stripe issue has not 
been specifically addressed in this paper. Through the 
microscopic separation of hole-rich regions from antifer- 
romagnetically correlated regions, |]6^ the stripe struc- 
ture tends to maintain the AF correlations more eff'ec- 
tively than the other case in which they arc uniformly 
suppressed by doped holes. Then it is not difhcult to 
expect that in the stripe state Tc and T* are suppressed 
by some amount from those values in the uniform state. 
We believe that the differences in the stripe state and 
in the uniform state are mainly quantitative. The con- 
clusion reached in this paper is not expected to qualita- 
tively change in the presence of stripe structure. How- 



ever, stripe-type structure can be important in the low 
doping region to stabilize both pairing and AF correla- 
tions at the same time. 

In some of magnetic experiments such as NMR, |6|] 
another crossover temperature r° (larger than T*) is of- 
ten identified, at which Knight shift shows its maximum. 
This feature can be easily understood on the basis of 
the competing nature of pairing and AF correlations as 
well as of the phase diagram obtained in Ref. The 
obtained mean-field phase diagram shows that the mean- 
field AF phase line stays above the mean-field SC phase 
line for x < Xc- As was discussed in Section HI, these 



AF and SC mean-field phase lines should be interpreted 
as the onset (T°) of short-range AF correlations and as 
the onset (T*) of short-range pairing correlations (pseu- 
dogap), respectively. For T* < T < T°, the two corre- 
lations compete, while they grow with decreasing tem- 
perature. In spin-lattice relaxation rate which picks up 
strongly the q = Q component, for T* <T <T^ the con- 
tribution from the AF correlations dominates that from 
the pairing correlations so that l/TiT keeps increasing 
until it starts to decrease at T* . On the other hand. 
Knight shift which picks up the (f = component and 
thus is unaware of the growing AF correlations, is more 
strongly influenced by the increasing pairing correlations. 
Consequently Knight shift reaches its maximum at T^, 
and starts to slowly decrease below it and then rapidly 
decrease below T* . 

Based on experimental results and our previous theo- 
retical work, we have conjectured a microscopic theory of 
high temperature superconductivity. In this conjecture, 
SC and AF long-range orders are driven by interlayer 
coupling. Strictly in two dimensions, the microscopic 
Hubbard model has an RVB insulator-to-metal transi- 
tion &i X = Xc near optimal doping for zero temperature, 
leading to a QCP, and one of the crossover lines is given 
by the pseudogap temperature T* . We argued that var- 
ious singular and non-Fermi liquid properties observed 
near optimal doping are due to the presence of this QCP. 
In our conjecture, the crossover line T* also practically 
divides the SC region into two, depending on the doping 
level with respect to Xc- For x < Xc the SC state has sig- 
niflcant AF correlations, while for x > Xc '^t has virtually 
no AF correlations, thus justifying the conventional BCS 
theory based on the noninteracting electrons. Inelastic 
neutron scattering resonance and systematically reduced 
superfluid density in the SC state below Xc have their nat- 
ural explanations in the present scenario. The present ap- 
proach supports interlayer pair tunneling model in which 
the SC condensation energy comes from the lowering of 
the c-axis kinetic energy in the SC state. Comparison 
of the present scenario with some of the leading theo- 
ries based on the Hubbard and i — J models was given. 
The generic features of both hole-doped and electron- 
doped cuprates as well as heavy-fermion superconductors 
may be understood in the unified framework within the 
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present picture. 
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FIG. 1. Schematic phase diagrams in doping (x = 1 — n) 
and temperature (T) plane based on (a) experiments and (b) 
our previous work, in the presence of the interlayer coupling. 
Tc and Tjv denote the SC and AF transition temperatures, 
respectively. T* is the pseudogap temperature, while Trvb 
in (b) is a temperature where the AF correlations become 
strong. The thick solid line in (b) denotes the insulating RVB 
ground state with a spingap. 
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FIG. 2. Schematic phase diagrams in doping {x = 1 — n) 
and temperature (T) plane (a) with a reduced interlayer cou- 
pling, and (b) and (c) without it. Tc, Tn, and T* are the same 
as in Fig. hi The thick solid lines denote the insulating RVB 
ground state with a spingap A, and without the interlayer 
coupling these lines terminate at x = Xc for T — shown as 
a black dot in (b) and (c). 
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